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Abstract 

We study the azimuthal asymmetry (cos 2(f)) in the Drell-Yan lepton pair production in hadronic 
scattering processes at moderate transverse momentum region, taking into account the contribu- 
tions from the twist-three quark-gluon correlations from the unpolarized hadrons. The contribu- 
tions are found to dominate the asymmetry, and are not power suppressed by q±/Q at small q± 
where q± and Q are the transverse momentum and invariant mass of the lepton pair. Accordingly, 
the Lam- Tung relation will be violated at this momentum region, and its violation depends on the 
twist-three functions. However, at large transverse momentum q± ~ Q, the Lam- Tung relation 
still holds because all corrections are power suppressed by A 2 /^ ~ A 2 /Q 2 where A is the typical 
nonperturbative scale. 
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1. Introduction. Drell-Yan lepton pair production in hadronic scattering process Jl[ has 



been playing a very important role in studying nucleon structure and QCD dynamics 



[2] 



and 



i i 

is an important complementary to the deep inelastic scattering studies [3J. Moreover, the 



later development on the angu 



and QCD dynamics studies 



ar distribution of lepton pair has laid ground for parton model 



6|. The lepton pair production in hadronic scattering, 



#1 + ij 2 -> 7* + x -> + x , (i) 



comes from the virtual photon decays. At higher energies, we should also consider the 
weak boson (Z°) decay contributions. In this paper, we will limit our discussions only 
for the virtual photon decays. An extension to including Z° boson decay contributions is 
straightforward. In the leading order, virtual photon isproduced through quark- ant iquark 



annihilation process, qq — * 7* in the parton picture [1]. In the rest frame of the lepton 
pair, we can define two angles [4]: one is the polar angle 9 between one lepton momentum 
and the hadron; the azimuthal angle is defined as the angle between the hadronic plane 
and the lepton plane. Here and in the following discussions, we follow the Collins-Soper 
frame {4] to define these angles. Our results can be translated to other frames too. The 



general expression for the lepton pair angular distribution can be written as 



^ = (1 + cos 2 9) + A (--- cos 2 9) + A 1 sin 29 cos + — sin 2 9 cos 20 . (2) 

It has been argued [4J that the coefficients A , A 1 , and A 2 are all power suppressed at 
large Q 2 , where Q is the invariant mass of the lepton pair: Aq ~ A2 ~ (k±) [Q 2 and 



k 



As a 



A\ ~ (k±)/Q where k± is the typical transverse momentum scale in the process 
result, the lepton pair angular distribution will be dominated by (1 + cos 2 9) in the small 
transverse momentum region. These power counting results were generalized to analyze the 
various relations between the above coefficients jlj]. One of the interesting observations is 
the so-called Lam- Tung relation Q]: 2v — (1 — A) = 0, where A = (2 — 3A )/(2 + A Q ) and 
v = 2A2/ (2 + Aq). According to the above power counting results, this relation is obviously 
valid because A = 1 and v — at the leading power. 

Finite transverse momentum of the lepton pair (q±) can be generated from gluon radiation 
from the leading partonic process, for example, through the quark-antiquark annihilation 
channel qq — > 7*g. Its contribution to the lepton pair angular distribution at finite transverse 



2 



momentum was found Q], 

dN 3 rg 2 + M q 2 -\qI 



+ 



COS 2 | 



sin 2 cos 20 + 



(3) 



^ i6vr L g 2 + gi g 2 + gi ~" " 1 g 2 

where the sin 26 term does not have simple expression and has been omitted in the above 
equation. Similar expression can be derived for the qg — > 7*g channel. Certainly, at finite 
transverse momentum, the angular coefficients Ai in Eq. (2) are not zero any more. However, 
they do obey the power counting rule. For example, the A 2 coefficient in Eq. (2) from the 
contribution in Eq. (3) is power suppressed by q]_/Q 2 at low transverse momentum. However, 
the Lam- Tung relation is still valid for any value of q± from this contribution. This has been 

higher order perturbative corrections to 
7], where it was found that the violation 



regarded clS cl simple prediction from QCD 
this relation has been calculated in the literature 



is numerically very small. These studies have motivated many experimental investigations 



of the lepton pair angular distributions in hadronic scattering 



In particular, the it 



induced fixed target experiment found large cos 20 azimuthal asymmetry, which is difficult 
to understand (J EJ. The intrinsic transverse momentum effects have been proposed to 
explain these effects [lo|], which however is limited to small transverse momentum whereas 
the experimental data are in both small and moderate transverse momentum region. 



Meanwhile, at low transverse momentum, there exist large 
a™ln 2n (g 2 /g 2 ) from the fixed order perturbative calculations 



11 



ogarithms in terms of 



12 



13j, Q Q. Re- 



summation of these larger logs has been formulated for the leading contribution from the 
partonic contribution, especially for the Drell-Yan lepton pair and vector boson production 



in hadronic scattering 14j]. This corresponds to the leading term in the lepton pair angular 
distribution (1 + cos 2 9) in Eq. (2). Because the rest terms (Ai) are power suppressed in 
the limit q±/Q -C 1, it has been difficult to follow the Collins-Soper-Sterman resummation 
method 9|. Recently, there have been much efforts to study the soft gluon resummation for 
these higher order terms [3] , and hopefully these developments will lead to a final solution 
to this issue, especially following the original Collins-Soper-Sterman formalism. 

In this paper, we study the lepton pair azimuthal asymmetry, in particular for the A 2 
coefficient in Eq. (2) from different perspective. We are interested in its behavior at the 
moderate transverse momentum region Aqcd € ?i < At this region, there are two 
large momentum scales q± and Q. The contribution to A 2 from Eq. (3) is power suppressed 
by q±/Q 2 as we mentioned above. However, from the following calculations we find that 



there exit contributions from the twist-three quark-gluon correlation functions from both 
incident hadrons, which are not power suppressed by qj_/Q 2 , instead by A 2 /qj_ where A is 
the typical nonperturbative scale. These contribution will dominate A 2 coefficient at the 
moderate transverse momentum region, depending on the relative strength of qj_/Q 2 and 
A 2 jq\ x . More importantly, because they are not power suppressed by q^/Q 2 , the soft gluon 
resummation can be performed following the classical Collins-Soper-Sterman approach, and 
the resummation effect will be similar to the leading term of (1 + cos 2 9) (see for example, 
the similar study in 17]). We notice that another higher-twist effect from only one side of 

nn 

the incoming hadrons has also been studied in the literature [18|, H9(, which are different 
from our calculations below. 

2. Twist-three times twist-three contributions to the lepton pair azimuthal 
asymmetry. From the general analysis of twist-three functions of the unpolarized 



hadrons 



20|, |2jJ, we find the on 



y twist-three function is T < f\x } y) which is equivalent to 



221 ]. It is defined as 



E(x,y) studied in the literature [2lL 

^i^2_^-( 3;2 -x 1 )P + +^r-i^ + (p|^(o-) ( T+^F + M (^)^(l/r)l^) 



F 



T£ , {x 1 ,x 2 ) 



47T 



(4) 



where \i is a transverse index, the sums over color and spin indices are implicit, \P) denotes 
the unpolarized hadron state with momentum P = (P + ,CT,0_|_) and P^ = (P° ± P z )/\/2, 
ip is the quark field, and F +fl the gluon field tensor, and the gauge link has been suppressed. 
This correlation is a chiral-odd function, and can generate transverse polarized Hyperon 
production in unpolarized hadronic collisions 23(|. Because of this chirality property, we 
have to introduce two correlation functions from both incoming hadrons. 
To calculate its contribution, we follow the 



velopments for the similar calculations 



25 



procedure outlined in 



26. 



20 



24] . and recent de- 



271 ] . In the collinear factorization framework, 



a general factorization formula for the contributions from the above correlation functions 
can be written as, 



da 



d^qdVt 



dxdx' dzdz' 



x 



(5) 



1 From power counting point of view, at this particular order, the new contribution will dominate in 
the region of A 2 <C q\ -C AQ whereas the contribution from Eq. (3) will dominate in the region of 
AQ « q\ « Q 2 . 
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where Tp(x,x') is the correlation function associated with the quark from hadron Hi and 



18 



20). 



T^(z, z') for the antiquark from hadron H 2} H. is the hard part and can be calculated from 
perturbative partonic process. This factorization formula follows (as a conjecture) earlier 
general arguments for the higher-twist contributions to the hadronic cross sections 
It will be very important to have a rigorous proof for this particular contribution as written in 
Eq. (5): the higher- twist effects coming from both sides of incoming hadrons, which is beyond 
the situations considered in [ijj]. Because of the higher-twist nature, it is always much 
more involved to calculate their contributions than those for the leading-twist contributions 
like Eq. (3). However, recent developments 25, 2^3, 27] have laid solid ground and useful 
technique to carry out those calculations. In particular, the two variables in will be 
fixed by taking pole contributions, or equivalently by calculating the imaginary part of the 
interference of the scattering amplitudes. For example, in the above equation x' will be 



equal to x depending on a soft or har e 
follow the procedure developed in 



25 



DO 



26 



e contribution 



25 



2fl 



271 ] . In this paper, we will 



271 ] to calculate the hard part in Eq. (5). 



First, we notice that the hard partonic part is separately gauge invariant summing up all 
possible diagrams. Therefore, we can carry out the calculations of these contributions with 
either A + or A± field connecting the hard and soft parts. A particular example has been 
given in (26!] for similar calculations. In our calculations, we find that it is more convenient 
to work the Aj_ part, and construct the field tensor accordingly. However, it has been known 
that the individual diagrams associated with the A± field depend on the boundary condition, 
although the final results of all diagrams contributions do not 3, 30]. Further calculations 
show that the retarded boundary condition A±(y~ = — oo) = will greatly simplify the 
derivations of Eq. (5). In the following, we will choose this boundary condition for the 
gauge field Aj_ from both hadrons. Under this boundary condition, on the other hand, we 
have to take into account the contributions from the operator with partial derivative on the 
quark field Njd±tp\ . This is because, this matrix element can be related to the quark-gluon 
correlation function defined in Eq. (5), 

dy 



ixP + y 



An 



Tp(.r) I (,:.r) 



F 



(6) 



with retarded boundary condition for the gauge field 29(. Therefore, for a complete calcu- 



lations, we have to ca 
with that of (ipA^ip) 



culate the diagrams associated with the operators (ipd±ip) , together 
2o| . For the A± contribution, we can further deduce its contribution 
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(a) 





FIG. 1: Generic Feynman diagrams for the twist-three times twist-three contributions to the lepton 
pair azimuthal asymmetry in hadronic scattering process: (a) d± contribution from both sides of 
hadrons: (b-c) d± and A± from each side; (d) A± from both sides. 

to that of Tp\xi.X2). For example, under the same boundary condition for the gauge field 
A i , we can write 



dy dy 1 
Air 



x 



P+ e - lP+ ^e t ^-^) p+ ^(P 1 S|^(0-)a + ^^ t (l/r)^(l/")l^) 
~Tp (x, X\) , 



(7) 



X\ + it 

where the pole structure in the second line comes from the partial integral and the ie 
prescription depends on the retarded boundary condition we are using. If we choose different 



y 



25 



29f . We have also checked 



26 



27] . 



boundary condition, this prescription shall change according 
that the above procedure can reproduce all previous results 

Following the above arguments, we plot the generic Feynman diagrams contributions in 
Fig. 1, where (a) is the contribution from d± operators from both sides of hadrons; (b-c) 
are those diagrams with d± and A± on either side; (d) for A± from both sides. There are 
four diagrams for Fig. 1(a), 12 diagrams for Fig. 1(b) and (c) respectively, and 208 diagrams 
for Fig. 1(d). As we mentioned above, in this paper we are interested in the cross section 
contributions in the moderate transverse momentum region q± <C Q. In carrying out these 
calculations, we will utilize the power counting method. We will only keep the leading power 
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contributions and neglect all higher power corrections in terms of q±/Q. The full expressions 
of our results will be presented in a separate publication. The advantage to use the retarded 
boundary condition is that we find that at the leading power, the contributions from the 
diagrams of Fig. 1(d) are either power suppressed or canceled out between soft and hard 
poles 29(| . We are left with the contributions from Fig. l(a-c) only, which are relatively 
easier to work out. In the limit of q± -C Q, the final result is 



da 



d A qd£l 



a sin 2 9 cos(20) 



gl2vr 2 



dx dz 

x z 



— — [ AT pi 



(z,z)5(C-l) + AT^(x,x)5(C-l) 



+2C F 5^-l)5(i-l)T^(x,x)T^-(z,z)\n% 

1± 



(8) 



where <r = a 2 m /6SQ 2 , £ = x /x, £ = z /z with x = -^e y and z = -S=e _y , and y is the 
rapidity of the lepton pair in the center of mass frame of incoming two hadrons and 5* is the 
hadronic center of mass energy square. The coefficient A is defined as 



A 



2N r 



X dx F « 



x, x) 



2Z + T^ q (x,x) 



2^-2) 



+ 



C A T {a), 



X,X Q ) 



(9) 



;i-0+ J ' 2 — 

and similar expression holds for A in the above equation. These diagrams (Fig. 1) also 
contribute to other terms including Aq and A\ in Eq. (2), but they are all power suppressed 
by qj_/Q 2 - From the above equations, we can see that this contribution to A 2 coefficient is 
not power suppressed by q\jQ 2 at moderate transverse momentum region. It is the same 
order as the (1 + cos 2 9) term in this power counting. Of course, it is suppressed by A 2 /q'j_ 
because of the higher-twist nature. This can also be seen from the above expression. 



More importantly, the above result can a 



dependent (TMD) factorization formalism 



Mulders function from both hadrons 
calculated in 



36 



28 



so be reproduced 



30. 



31 



32 



33 



34 



yy a transverse momentum 



35] with the so-called Boer- 



37] at large transverse momentum which has been 



231 ] . This demonstrates that in the intermediate transverse momentum region, 



for this part contribution, the twist-three times twist-three collinear factorization approach 
and the TMD factorization approach are consistent for the cos(20) azimuthal asymmetry 
in the unpolarized Drell-Yan processes. This is a nontrivial demonstration, because it goes 



beyond previous examples studied in the literature 27j where the twist-three effect from 
only one side of the incoming hadrons was considered. Because of this consistency, the 
energy evolution equation [13j (the Collins-Soper evolution equation) can be derived for 



this contribution, and the soft gluon resummation can be accordingly performed. This will 
significantly change the relative sizes of this contribution and the contribution from Eq. (3) 



to the cos 20 asymmetry. We will leave a detailed study in a separate publication 29]. 

3. Conclusion. We have the following results for the angular distribution of the Drell- 
Yan lepton pair production in hadronic reactions, 

• At moderate transverse momentum, A 2 is in order of 1, A Q is power suppressed by 
q\/Q 2 . As a result, the Lam- Tung relation will be violated because A is 1 whereas v is 
order of 1. Of course this violation will depend on the sizes of the twist-three correlation 
function from both incoming hadrons. Furthermore, soft gluon resummation 
will not change the power counting result for A 2 , because the leading order TMD 
factorization leads to the same resummation pattern similar to that discussed in jl4|. 
which is important to understand the angular distributions of the lepton pair at this 



momentum region 



9]. 



• At large transverse momentum q± ~ Q, however, both Aq and A 2 are in order of 1, 
and they are dominated by the contribution from the unpolarized quark and antiquark 
annihilation contribution Eq. (3). The contributions we calculated this paper are 
suppressed by A 2 jq\ ~ A 2 /Q 2 . Because of this, the Lam- Tung relation will be valid 
again. 

In summary, we have investigated the higher-twist effects to the Drell-Yan lepton pair 
angular distributions in hadron-hadron scattering processes. We found that the twist-three 
times twist-three contributions to the cos 20 azimuthal asymmetry contribution are not 
power suppressed by q^/Q 2 , rather by A 2 jq\ at the moderate transverse momentum. We 
further argue that this part of contribution will not be affected by the soft gluon resumma- 
tion effects, and the Lam- Tung relation will be modified at small and moderate transverse 
momentum. It will be interested to compare our predictions with the experimental data 8J 
and check the phenomenological importance of our results. It will also be interested to 
extend to other processes like the cos 20 asymmetry in semi-inclusive hadron production 
in deep inelastic scattering and back-to-back two hadron production in e + e~ annihilation 
processes where the similar effects shall play very important roles. 

We thank Jianwei Qiu and Werner Vogelsang for interesting discussions. This work was 
supported in part by the U.S. Department of Energy under contract DE-AC02-05CH11231 

8 



and the National Natural Science Foundation of China under the approval No. 10525523. 
We are grateful to RIKEN, Brookhaven National Laboratory and the U.S. Department of 
Energy (contract number DE-AC02-98CH10886) for providing the facilities essential for the 
completion of this work. J.Z. is partially supported by China Scholarship Council. 



[1] S. D. Drell and T. M. Yan, Phys. Rev. Lett. 25, 316 (1970) [Erratum-ibid. 25, 902 (1970)]. 
[2] R. Brock et al. [CLEQ Collaboration], Rev. Mod. Phys. 67, 157 (1995). 
[3] E. D. Bloom et al, Phys. Rev. Lett. 23, 930 (1969); M. Breidenbach et al, Phys. Rev. Lett. 
23, 935 (1969). 

[4] J. C. Collins and D. E. Soper, Phys. Rev. D 16, 2219 (1977). 

[5] C. S. Lam and W. K. Tung, Phys. Rev. D 18, 2447 (1978); Phys. Lett. B 80, 228 (1979); 

Phys. Rev. D 21, 2712 (1980). 
[6] J. C. Collins, Phys. Rev. Lett. 42, 291 (1979). 
[7] E. Mirkes and J. Ohnemus, Phys. Rev. D 51, 4891 (1995). 

[8] See for example, J. S. Conway et al, Phys. Rev. D 39, 92 (1989); L. Y. Zhu et al [FNAL- 

E866/NuSea Collaboration], Phys. Rev. Lett. 99, 082301 ( 2007k laFXTv:Q811.4589 l [nucl-ex]. 
[9] P. Chiappetta and M. Le Bellac, Z. Phys. C 32, 521 (1986). 
[10] D. Boer, Phys. Rev. D 60, 014012 (1999). 

[11] Y. L. Dokshitzer, D. Diakonov and S. I. Troian, Phys. Lett. B 78, 290 (1978); Phys. Rept. 
58, 269 (1980). 

[12] G. Parisi and R. Petronzio, Nucl. Phys. B 154, 427 (1979). 

[13] J. C. Collins and D. E. Soper, Nucl. Phys. B 194, 445 (1982). 

[14] J. C. Collins, D. E. Soper and G. Sterman, Nucl. Phys. B 250, 199 (1985). 

[15] D. Boer and W. Vogelsang, Phys. Rev. D 74, 014004 (2006). 

[16] E. L. Berger, J. W. Qiu and R. A. Rodriguez-Pedraza, Phys. Lett. B 656, 74 (2007); Phys. 

Rev. D 76, 074006 (2007). 
[17] A. Idilbi, X. d. Ji, J. P. Ma and F. Yuan, Phys. Rev. D 70, 074021 (2004). 
[18] J. W. Qiu and G. Sterman, Nucl. Phys. B 353, 105 (1991); Nucl. Phys. B 353, 137 (1991). 
[19] V. Braun, R. J. Fries, N. Mahnke and E. Stein, Nucl. Phys. B 589, 381 (2000) [Erratum-ibid. 

B 607, 433 (2001)]. 



9 



[20] R. K. Ellis, W. Furmanski and R. Petronzio, Nucl. Phys. B 212, 29 (1983). R. K. Ellis, 

W. Furmanski and R. Petronzio, Nucl. Phys. B 207, 1 (1982). 
[21] R. L. Jaffe and X. D. Ji, Nucl. Phys. B 375, 527 (1992). 
[22] J. Zhou, F. Yuan and Z. T. Liang. larXiv:0812 .44841 [hep-ph]. 
[23] J. Zhou, F. Yuan and Z. T. Liang, Phys. Rev. D 78, 114008 (2008). 

[24] J.W. Qiu and G. Sterman, Phys. Rev. Lett. 67, 2264 (1991); Nucl. Phys. B 378, 52 (1992); 

Phys. Rev. D 59, 014004 (1998). 
[25] C. Kouvaris, J. W. Qiu, W. Vogelsang and F. Yuan, Phys. Rev. D 74, 114013 (2006). 
[26] H. Eguchi, Y. Koike and K. Tanaka, Nucl. Phys. B 752, 1 (2006); Nucl. Phys. B 763, 198 

(2007). 

[27] X. Ji, J. W. Qiu, W. Vogelsang and F. Yuan, Phys. Rev. Lett. 97, 082002 (2006); Phys. Rev. 

D 73, 094017 (2006); Phys. Lett. B 638, 178 (2006). 
[28] X. Ji and F. Yuan, Phys. Lett. B 543, 66 (2002); A. V. Belitsky, X. Ji and F. Yuan, Nucl. 

Phys. B 656, 165 (2003). 
[29] J. Zhou, F. Yuan, and Z. T. Liang, to be submitted. 

[30] S. J. Brodsky, D. S. Hwang and I. Schmidt, Phys. Lett. B 530, 99 (2002); Nucl. Phys. B 642, 
344 (2002). 

[31] J. C. Collins, Phys. Lett. B 536, 43 (2002). 

[32] D. Boer, P. J. Mulders and F. Pijlman, Nucl. Phys. B 667, 201 (2003). 

[33] X. Ji, J. P. Ma and F. Yuan, Phys. Rev. D 71, 034005 (2005); Phys. Lett. B 597, 299 (2004). 
[34] J. C. Collins and A. Metz, Phys. Rev. Lett. 93, 252001 (2004). 

[35] P. J. Mulders and R. D. Tangerman, Nucl. Phys. B 461, 197 (1996) [Erratum-ibid. B 484, 
538 (1997)]. 

[36] D. Boer and P. J. Mulders, Phys. Rev. D 57, 5780 (1998); D. Boer, S. J. Brodsky and 

D. S. Hwang, Phys. Rev. D 67, 054003 (2003). 
[37] A. Metz, P. Schweitzer and T. Teckentrup. larXiv:0810^2T2l [hep-ph]. 



10 



